Multiple light scattering in nematic liquid crystals 
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We present a rigorous treatment of the diffusion approximation for multiple light scattering in 
anisotropic random media, and apply it to director fluctuations in a nematic liquid crystal. For a 
typical nematic material, 5CB, we give numerical values of the diffusion constants and D±. We 
also calculate the temporal autocorrelation function measured in Diffusing Wave Spectroscopy. 
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Light transport in random or turbid media has long 
been treated by radiative transfer theories, the first of 
which was formulated as early as 1905 by Schuster 0. 
For distances large compared to the transport mean-free- 
path I*, beyond which the direction of light propagation 
is randomized, these theories can be reduced j|] to a dif- 
fusion equation for the light energy density with diffu- 
sion constant D = cl* /3 where c is the speed of light in 
the medium. In 1984 Kuga and Ishimaru Q discovered 
coherent backscattering of light in colloidal suspension, 
predicted in earlier papers Q| , and physicists realized the 
connection of wave propagation in disordered media to 
weak localization jsj , a precursor of Anderson localization 
Since then, our theoretical understanding of light 
transport in random media has advanced considerably. 
Detailed studies of multiple scattering of scalar waves |^] 
was followed by the generalization to include the polariza- 



tion of light 
cal activity 




, broken time reversal symmetry and opti- 
and long-range correlations in random 
scatterers |l^ ll|] . Multiple scattering emerged as a pow- 
erful probe of dynamical properties of turbid media with 
the development of Diffusing Wave Spectroscopy (DWS) 
[ p2| , p^ as an experimental technique capable of measur- 
ing dynamic correlations at time scales, much shorter 
[Q , than can be probed with single scattering. 

Nematic liquid crystals are strong light scatterers, ex- 
hibiting turbidity and coherent backscatter [p^ . They 
differ, however, in significant ways from colloidal suspen- 
sions, the most widely studied multiple-scattering media. 
First, nematic liquid crystals are anisotropic with bar like 
molecules aligned on average along a unit vector n{r,t) 
called the director. They are birefringent with different 
velocities of light for ordinary and extraordinary rays. As 
a result the photon energy density, like particle density 
in an electron system [p6| , obeys an anisotropic diffu- 
sion equation with diffusion coefficients and D± for 
directions parallel and perpendicular to the equilibrium 
director Uq. Second, the dominant scattering of visible 
light is from long-range thermal fluctuations of the direc- 
tor rather than from particles with diameters comparable 
to the wavelength of light. This leads to a divergent scat- 
tering mean-free-path when the external magnetic field 
H is zero |23 2^. The diffusion constants D|| and D±, 



which in isotropic systems are proportional to the trans- 
port mean-free-path, are, however, finite when H 
as shown in Fig. 0. In this paper, we develop a system- 
atic treatment of the diffusion approximation for multiple 
light scattering in anisotropic random media (for recent 
approaches see jl^), which allows us to calculate D^^ and 
D± from known properties of nematics and to obtain the 
time-dependent response measured in DWS experiments 
|]l8| . Fig. |l| shows our calculated values of Z3|| and D± 
as a function of external magnetic field H for the com- 
pound 5CB. The anisotropy ratio D\\/D± — 1.45 is in 
good agreement with measurements reported in a com- 
panion paper [19] . 
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FIG. 1. The field dependence of the normalized diffusion 
constants Dy and D± and the anisotropy (_D|| — D±)/D± 
for parameters of a typical nematic liquid crystal 5CB: 
K1/K3 = 0.79, K2/K3 = 0.43 and Ae/e± = 0.228. 

We start with the wave equation for the electric light 
field E{r,t): 
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[£o + 6eir,t)]'}E{r,t)^0 . (1) 



The homogeneous part of the dielectric tensor is £o; Sind 
the randomly fluctuating part ^£(r, t) is a Gaussian ran- 
dom variable described by the correlation function 



B'^{R,t) 



[Se{R,t))(g)Se{O,0) 



(2) 
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where lo is the frequency of hght. The superscript (iV) 
means that we interchange the second and third index 
in the tensor product 5e ® Se to define B"^ : [B'^]ijki cx 
( Ssik Ssji ). We cah B^{R, t) the structure factor of the 
system. It is measured in single light scattering exper- 
iments p^ l and contains informations about the elastic 
and dynamic properties of a system. The local uniaxial 
dielectric tensor can be expressed as 



£(r, t) = e±_l + Ae[n(r, t) (g) n{r, t)] 



(3) 



Here e± and ey are the dielectric constants for electric 
fields, respectively, perpendicular and parallel to the di- 
rector, and Ae — £|| — e±. We assume that the inho- 
mogeneity of the director field only comes from thermal 
fluctuations of the director around its equilibrium value 
no: n{r,t) = tiq + 6n{r,t), where Sn has to be per- 
pendicular to no for small fluctuations. The dominant 
contribution to {R, t) is proportional to the director 
correlation function ( 5n{R^ t) (g) dn{0, 0) ) which we ex- 
press in momentum space |2l| ]: 



[6n{q,t)(g6n*{q,0)) = 



Ek^T 



KM) 



Ua{q) ® Ua{q) ■ (4) 



Here Ka{q) = I<aqj_+K3q'^^+AxH^, where Ki, K2, and 
are the Frank elastic constants, H is the external field 
parallel to no, and Ax is the anisotropy of the magnetic 
susceptibility. The quantity rja{q) is a combination of 
viscosities which appear in the hydrodynamic equations 
of the director field, the Leslie-Erickson equations 
The unit vectors Ua{q) specify the direction of Sn{q,t) 
in mode a — 1 and 2. 

In an anisotropic medium with a homogeneous dielec- 
tric tensor £0 the electromagnetic field travelling along 
the unit vector k has two modes with indices of refrac- 
tion na{k) and electric polarization ea(k). The polariza- 
tion tZ"(fc) of the displacement field, d"(fc) = £0^0 (fe), 
obeys d°'{k) ■ k — 0. After an appropriate normal- 
ization, the vectors fulfill the biorthogonality relation 
d"{k) ■ ef}{k) — dp [Q. Wc can now write down the mo- 
mentum space representation of the averaged retarded 
and advanced Green's function of Eq. in the weak- 
scattering approximation: 



(fc)} 'e„(fc)®e„(fe) (5) 



with 



la{k,Uj) 



(6) 



being the scattering mean-free-path of a light mode 
{k°' \ ea{k)} with wave vector fc" = ^Uak and polar- 
ization ea{k). The structure factor = 0)]q/3 

describes the scattering of {fe" | ea{k)} into {q'^ \ ep{k)}. 
The symbol J.^ always stands for an angular integration: 



(27r) 



(g) 



(7) 



To the order of our calculations ( G^^^ ) (fc, lo) is diagonal 
in the polarization ea{k). In what follows, a Greek index 
will refer to the "basis vectors" Cq ® Ba or d" g)d°'. The 
scattering mean-free-path la{k,u}) in the nematic phase 
has been calculated (2^j2^ . For small H for the extraor- 
dinary ray it has the form l'^ cx (|)'^ln(^^), 
where K is an appropriately averaged elastic constant. 

Let us look at the spatial and temporal autocorrela- 
tion function for the electric light field: ( E{R + + 
^) ® E*{R— — |)), where we have already intro- 
duced center of "mass" {R,T) and relative {r,t) co- 
ordinates. From this quantity, others follow as spe- 
cial cases: the energy density of light at time T is 
Wi{R,T) = {E{R,T) ■ eQE*{R,T)), where the T de- 
pendence is, e.g., due to time dependent-sources; the 
temporal correlation function of a steady-state light field 
is W2{R,t) = {E{R,\) ■ £oE*{R,-^)), which refiects 
the dynamics of the scattering media measured in DWS 
experiments. The Fourier transform with respect to r 
gives the energy density with wave vector A; |^ . To calcu- 
late the autocorrelation function for special light sources 
and/or given boundary conditions, we need the "two- 
particle" Green's function * = (G^®G^)W. Our 
goal is to derive the diffusion pole of $ in momentum 
and frequency space. With all arguments, the Green's 
function is {K, 57, t). K, correspond to the center 
of "mass" coordinates R, T and k, k' to the relative co- 
ordinates r,r' . The superscript lo is the light frequency, 
and the t dependence explicitly comes from the struc- 
ture factor B^j., (i). In the weak-scattering approxima- 
tion, (JC, £7, t) can be represented as a sum of lad- 
der diagrams, which is equivalent to the Bethe-Salpeter 
equation: 



/ ^[i^^l-mK,n)B%^^{m%^^,{K,n,t) 



where 



mK,n)^[{G"){k+,u;+)® (G-^)(fc_,c^_; 
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(8) 



(9) 



with k± = k±K/2 and a;± —Lu±fl/2 and, [i-l^lr]ijki ■= 
{2tt)^ S{k — k') {SikSji + SuSjk)/2. The multiple integrals 
and the sum can be done analytically for (5-function cor- 
relations but not for the anisotropic, long-range correla- 
tions of our problem. VoUhardt and Wolfle pst] derived 
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the diffusion pole for isotropic electron transport directly 
from Eq. (||), and MacKintosh and John applied their 
method to light. In the anisotropic case one has to be 
more careful |l6|. If ^^^\K,n,t) and X^''\K,n,t) are, 
respectively, the nth eigenvector and eigenvalue of the 
integral operator of Eq. (||), 



d^ki . (4) 
(27r)3 ^ 



(») 

k 



(10) 



and {K, 57, t) the eigenvectors of the Hcrmitian ad- 
joint operator, it is straightforward to show that 

{K,n,t) = Y. ^^nf"' n'iK,n) (11) 

n 

solves the Bethe-Salpeter equation [0. In the case of 
K = and Q = t = 0, it can be shown that the quantity 
AG^(O,0), where 

AG^(K,fi) = (G«)(fe+,c.+)-(G^)(fc_,u;_) , (12) 

is an eigenvector with eigenvalue A*^°^ (0,0,0) = 0. This 
is a very general result, based on the Ward identities 
valid beyond the weak-scattering approximation [ p5| . We 
have identified the diffusion pole, as we shall explicitly 
see soon. All other eigenvalues are positive, and, in real 
space, they give exponentially decaying contributions to 
4'^^,(fC, 17, i) of Eq. ( pT| ) , which are not important at 
long length scales |jll[ . To establish the diffusion approxi- 
mation we have to apply perturbation theory to calculate 
^(0) ( f'j fQj. sniall K, il, and t. Therefore, we expand 
the eigenvectors into a set of basis functions and turn the 
eigenvalue equation ( |lo|) into a matrix equation. For the 
component a of 'S'fc we use the ansatz 

vI/^(x[AG^(0,0)]"[$o7r + ^$?<(fc)] , (13) 



where 



LO \ C 



V c 



(14) 



The first factor on the right-hand side of Eq. (|13|) is due 
to the momentum shell approximation, Eq. (|l4[); it is 
strongly peaked around the wave numbers of the light 
modes. The amplitude ^'o then represents the zeroth 
eigenvector. The second term corresponds to the space of 
all other eigenvectors where the Lpf (k) are basis functions 
on the unit sphere, e.g. spherical harmonics, which can 
in general depend on polarization a. We also use the 
relation 



X AS^(O,0) 



[/fc(O,0)]° 
dGn 



dk 



-K 



(15) 



which gives [AG^(-K', i7)]" correctly to first order in 

and S (see [^). (S is the mass operator and AS^ 
is defined the same way as AG^. Gq stands for the 
Green's function of the homogeneous medium.) The cou- 
pling between the zeroth and the other eigenvectors then 
produces the diffusion tensor, and, finally, the Green's 
function takes the form: 



^tk'{K,n,t) 



1 ag^(o,o)0Ag;^,(o,o) 

N -in + fj.{u;, t) + K ■ D{uj)K 



(16) 



with N — —2n'^u}'^/{TTc^) and being the angular and 
arithmetic average of the two refractive indices. The de- 
nominator represents a diffusion pole, which also contains 
an "absorption" coefRcient ii{uj,t). The diffusion tensor 
follows from 



K ■ D{u;)K = -= [g{K)]* ■ B-'g{K) (17) 
2n-' 



with 



[g{K)U = vr / n„(fe) 



(fe) [v7{k)Y 



and 



^ ^ Jk" Jq-/ 

-Mik)r^^iq)s}]m.^-,{o)u] . 

In principal all ff{k) of odd parity contribute to D{lo). 
For isotropic systems we choose spherical harmonics: 

(fc) Yim{'d,'p)- Only the components [g{K)] 
are nonzero and [S]^;™^, oc Su'. Therefore, only spher- 
ical harmonics of ^ = 1 contribute to D{uj) and we get 
the familiar formular D — ^cl* oc [( 1 — cos'iS)]"^. The 
absorption coefRcient reads 



M(a;,t) = -= 2^ 



k''qP 



It represents an angular average over all the dynamical 
modes of the system. (For i = 0, it is zero and then 
increases due to the decaying temporal correlations in 
( (5e (Je ).) The numerator in Eq. (|l^) indicates which 
initial and final polarization states have a nonzero overlap 
with the diffusion pole. The second factor AG^/ (0, 0) de- 
pends only on the input wave. The first factor AG^(0, 0) 
involves only the output wave and determines the ra- 
tio of densities of photons in the two output polariza- 
tion states 1 and 2 independent of the state of the input 
wave. An integration over k (/ k'^dk) shows that this ra- 
tio is [ni{k) / n2{k)]^ for the wave direction k. This effect 
should be measurable. Finally, the Green's function cor- 
responding to W2{R,t) follows from ^'^^,{K,n = 0,t) 
by integrating over fc, k' and applying the appropriate 
trace operation. 
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The diffusion tensor D{uj) has the same uniaxial form 
as the dielectric tensor in Eq. (||). We express the diffu- 
sion coefficients Dh and D± in terms of a typical length 



/* 



97r' 



kuT As 



ical factors D 



and D i 



c/y^eX) times unitless numer- 



via 



(18) 



is over time r rather than path length because the light 
velocity is not a constant along an arbitrary path.) 
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where and D± depend on K1/K3, K2/K3 and Ae/e±. 
For the material 5CB, K3 = 5.3 x IQ-'^dyne, £± = 2.381 
and Ae/ejL = 0.228. With T = 300 K and green 



light (w/c = 1.15 X 10^ cm ^) we get IX 



2.3 mm, 



which is in agreement with experiments P,|l§. As ba- 
sis functions we choose spherical harmonics depending 
on a new "polar angle" -d'^ (see |^). For the ordi- 
nary light ray i?2 — i?, for the extraordinary one "i?^ 
is given by cos{>[ — ni{k)cos'd / y^WT- The basis func- 
tions ipfjy^{k) = Yimi^aii^) , are orthogonal with re- 
spect to the weight Til^{k)dflk oc dcos-d'^dip. Then, only 
[G{K)]ai=im is nonzero |2^. We studied the contribu- 
tions of different I to _D|| and D± and found that I = 3 
in addition to / = 1 gives changes of less than 2%. In 
Fig. ^ we plot our results for Z)|| and D± for 5CB with 
K1/K3 = 0.79, K2/K3 = 0.43 and Ax = 0.95 x lO"'^. 
At if = 0, Z?!! = 0.95 and D± = 0.65 are finite even 
though, as noted earlier, the scattering mean-free-path 
for the extraordinary light ray is infinite. The anisotropy 
in the diffusion constants decreases with both As and 
anisotropy in the Frank elastic constants. In the limit 
Ae = and Ki ^ K2 = K3, D\\/D^ = 1.06 is not unity 
because of the inherent anisotropy in the structure fac- 
tor. The diffusion approximation is valid only for times t 
much smaller than characteristic relaxation times of the 
director modes. In this case we get 



2fc, 



97r 



5T Ae2 jl 
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(19) 



where 7 is the rotational viscosity and jl a numerical fac- 
tor depending on all other viscosities and Ae/ej_. Note 
that unlike scattering in colloids, depends only on vis- 
cosities and is independent of the static structure factor 
(oc UbT / Kq^). This is because the same fluctuations 
determine scattering and dynamics. Finally, we point 
out that the appropriate Laplace-Fourier transform of 
Eq. (16) leads to a temporal autocorrelation function W2 
that can be expressed in a form reminiscent of the aver- 
age over light paths used in isotropic systems |l2^,f3| : 



IV2 oc / dr P{t) exp{~fiotT) 



(20) 



where P{t) is the probability that an anisotropic random 
walker enters the medium at a prescribed point and leaves 
it at another point after a time r. (Note that this integral 
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